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Abstract 





qh| The standard experimental techniques usually adopted in the study of the behaviour 

of ultracold atoms in optical lattices involve extracting the atom density profile from 
absorption images of the atomic sample after trap release. Quantum mechanically this 
I procedure is described by a generalized measure (POVM); interference patterns found 

in absorption images suggest a generalized measure based on fixed-phase, coherent-like 
(J\ . states. We show that this leads to an average atomic density which differs from the 

O I usually adopted one, obtained as the expectation value of the atom density operator 

in the many-body state. 

OO 
O 

> 

X ■ 1 Introduction 



Experiments with ultracold atoms trapped in optical lattices have reached high degrees of 
accuracy and sophistication in exploring the coherence effects that such systems exhibit; 
they turn out to be a unique laboratory for experimental tests of many effects in many-body 
quantum physics and in particular in the study of quantum phase transitions (for recent 
reviews see [H E] , and references therein) . 

In all these studies, the technique used to extract information on the system is essentially 
the same and is based on the analysis of interference phenomena (see [3]- [15] and references 
therein). Being rather difficult to measure relevant system observables directly inside the 
optical lattice, a commonly adopted experimental procedure consists in releasing the con- 
fining trap and subsequently obtaining absorption images of the expanding atomic clouds. 
Since during expansion the atoms are usually considered non interacting and thus following 
a free, ballistic evolution, the set of the obtained absorption images contain information on 
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the dynamics of the atoms in the lattice prior to the trap release. The imaging process is 
destructive: the sample is lost after the measurement. Then, one usually repeats the whole 
procedure, preparing the system in the same initial state, and obtaining at the end a set of 
absorption images. When all images are put together and superimposed, in the limit of a 
large number of them, they are usually assumed to reproduce the atom density obtained by 
averaging the associated operator in the given many body system state. 

In general, the quantum mechanical description of the averaging procedure corresponds to 
a generalized measure, a so-called Positive Operator Valued Measure (POVM) [161 [TTl [H] . 
each single image obtained by "taking a photograph" of the expanding clouds, the process 
of photon detection gives rise to interference patterns showing a phase difference among 
different wells, even though the initial state does not have a definite phase difference: the 
imaging process effectively projects the state of the system into a coherent-like state, with 
definite phase relations, although varying randomly from picture to picture [121 E]- Given 
this, one can compute using the corresponding POVM the atom density profile: we show 
it to be in general different from the one obtained through a direct average of the operator 
density. 

In most cases this difference turns out to be negligible due to the large number of atoms 
involved. Nevertheless, at least in principle, it can be made apparent by preparing the 
system in suitable states, for which the number of atoms in the various lattice sites is highly 
unbalanced. In this respect, a very promising instance in which this difference might be 
experimentally investigated is provided by the measure of density-density correlations in 
two-color optical lattices [191 120] • 

In the following, we shall discuss the case of an optical lattices formed by just two sites, 
i. e. a system of ultracold atoms confined in a double- well potential: besides simplifying the 
discussion, this case has direct experimental relevance; generalization of our considerations 
to a generic lattice presents no difficulties. 

We shall first briefly recall the two-mode description of the trapped system and then 
analyze its behaviour after the release of the trapping potential using a second-quantized 
many-body approach. The average density profile obtained by superimposing the various 
absorption images will then be computed starting from an initial Fock state, containing 
a definite number of atoms in each well, assuming, as often done in the literature, a free 
expansion of the atomic cloud after trap release. We shall see that for a very unbalanced 
initial state, this average differs in an experimentally testable way from the standard density 
mean value. 

2 Cold bosonic gas in a double well trap 

In a suitable approximation, i.e. for large enough barrier, the dynamics of cold atoms 
confined in a double well potential can be described by a two-mode Bose-Hubbard type 
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Hamiltonian 

H = El a\ai + E2 ala2 + U {a\Y a\ + (og)^ a\ — T {a\a2 + aia. 
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where 01,2, 2 annihilate and create atoms in the first, respectively second well and sat- 
isfy the Bose commutation relations [a^ , aj] = 5ij. The three contributions in the above 
hamiltonian correspond to a hopping term dependent on the tunneling amplitude T, and 
on two on-site energy terms, one quadratic in the number operator with coupling constant 
[/, describing the boson-boson repulsive contact interaction, the other due to the trapping 
potential, with the parameter E2 representing the wells depth. 

Notice that the total number of particles is conserved by the Hamiltonian ([1]). As a 
consequence, the Hilbert space of the system is + dimensional and can be spanned by the 
set of Fock states, describing the situation in which the first well is filled with k particles, 
while the other contains N — k atoms; they are generated by the action of the creation 
operators on the vacuum: 

Alternatively, one can introduce coherent-like states, depending on two real parameters, 
an amplitude ^ G [0, 1] and a phase ip G [0, 27i], 

\N;^,V) = (v^e^^/^al + v^e-^^/^a^)"" |0) . (3) 

They describe situations in which all A^ particles are in a coherent superposition, with definite 
relative phase <f and relative mean occupation number (A^; ^, ip\ (ajoi— 0302) |A^; ^, <f) = 1—2^. 
The set of states ([3]) is overcomplete, 

and although not strictly orthogonal, they become so in the limit of large A^: 

(AT; e, ifiN; e', ^ - o - V') ■ (5) 

Nevertheless, one can always express any Fock state ([2]) in terms of coherent states by using 
the completeness relation (jlj, 

I AT; k) = {N+ 1) r ^ {N; e, AT; k) \N; ^, ^) , (6) 

where the overlap functions are explicitly given by 



(AT; e, ^\N- k) = ( V ) (1-0^ e-^H'=-T) . (7) 
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As well known [T], the Hamiltonian ([T]) describes a cross-over between a superfluid and 
insulator phases, which becomes a true quantum phase transition, with order parameter 
depending on the ratio T/f/, in the limit of an infinite number of wells. Neglecting the shift 
E = El — E2, when an atom hops between the two wells, it loses an energy T and gains 
instead an energy U; thus, when T/U <^ 1, this situation is energetically suppressed and, at 
equilibrium, the ground state of the system is given by a Fock state with equal number N/2 
of bosons per well, |A^; (Mott insulator phase). On the other hand, when T/U ^ 1, 

the system shows phase coherence; all the N particles are in the same superposition and 
the ground state of ([1]) can be approximated by \N; ^, ip), with a definite relative phase and 
occupation number (superfiuid phase). 



3 Generalized quantum measure 

As explained in the introductory remarks, in a typical experimental set up the double well 
trap is filled with a large number of cold atoms, whose dynamics is driven by the Hamilto- 
nian ([1]). Since in-trap measures of relevant system observables are difficult, one can obtain 
indirect information on the atom behaviour by switching off the trapping potential and let- 
ting the atoms expand freely. The two fractions of atoms coming from the two wells will 
eventually overlap; one is then able to observe the spatial distribution of the atoms by il- 
luminating the overlapping clouds with a probe light and by registering the corresponding 
absorption image. 

In order to model this process of measure, we shall use a second-quantized many-body 
formalism. Let us first introduce a complete set of single-particle atom states 
obtained by acting with the set of all creation operators a] on the vacuum, \Wi) = al\0). 
Although just two creation operator a[ and a\ are enough to properly describe an atom 
confined into the two wells, in absence of the trapping potential a single atom moves freely 
in space and therefore a complete set of wave functions Wi{x) = {x\Wi) = {x\al\0) is needed 
for properly representing its states [x represents the collection of spatial coordinate vari- 
ables). The creation operator ip'^i^x) of an atom at position x, ip'^{x) |0) = |a;), can then be 
decomposed as 

00 

ij\x) = J2w,{x)al . (8) 

i=l 

Since the single particle states are orthonormal, one can invert this relation and equivalently 
write 

a] . j ,.W.(.m^) . ,9) 

Further, from [a^, a]] = = Sij, one recovers the standard bosonic (equal-time) 

commutation relations: 

[ip{x), ^p\x')] = 6{x-x') . (10) 
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With this formahsm, one can now easily compute the action of the operator iIj{x) on a 
many-body state |A^;^,<^). Using the above commutation relations, one easily finds: 

ij{x)\N-^,^) = Vn (^^/^W,{x)e'^/^ + ^/T^W^ix) e''^/^^ \N-1;^,^) , (11) 

so that the wave function associated to the coherent state \N;^,ip), depending on the posi- 
tions Xi, X2, ■ ■ ■ , xn of the atoms, is simply given by 

N 

{xi,X2, . . . ,xn\N-^,^) = Vn^. l[(^^Wi{x,)e"^/^ + y^T^W2{x,)e-"^/^^ . (12) 

1=1 

Similarly, the action of the destruction operator iIj{x) for an atom at position x on a Fock 
state |A^; k) can be obtained from (fTTj) through the expansion IQ. 

Strictly speaking, all this holds at fixed time, e.g. at t = 0. However, since the atoms are 
assumed to evolve freely once the trap is released, insertion of the time dependence in the 
formalism is straightforward. Assume that the confining double well potential is released at 
time t = and denote by Ut the unitary operator that evolves freely in time the initial one 
particle states: 

\W,{t)) = U^m = Utal\0) = al{t)\0) , al{t) = Uta\ul . (13) 

The corresponding wave function is given by Wi{x] t) = {x\Wi{t)) = (a;|a|(t)|0), and coincides 
with a transformed Wi{x) under a ballistic expansion. Since the dynamics is free, every 
particle in a many-body state will evolve independently with Uf, therefore, the evolution up 
to time t of the t = coherent state \N; ^, ip) will simply be given by 

\N;^,^;t) = (v^e^^/^aKt) + 6-^/^4 (t))"" |0) , (14) 

and similarly for a Fock state (I2l). 

In this picture, states evolve in time while observables remain fixed and thus the operator 
tp{x) results time-independent. Since at each instant of time t the collection {|VFi(t))}^i is a 
complete set of single particle states obtained from the vacuum by the action of the creation 
operators a|(t), ipi^) can be equivalently decomposed as ip{x) = VFj(x; t) aj(t) for all 

times. The action of ip{x) on \N;C,,<f',t) is then as in (fTTj) . once expressed in terms of the 
evolved wave functions: 

ij{x)\N-^,^-t) = Vn (^^W,{x-t)e'^/^ + ^/r^W2{x-t)e-''^/^^ \N-l-^,^-t) . (15) 
Through this result, one can easily compute the mean value of the density operator at point 

X, 

n{x) = ijj^x) ^p{x) , (16) 
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in any many-body state. For instance, if the atoms just before the release of the confining 
potential are in a superfluid state described by a state |A^; ^, ip), after a time t of free evolution 
of the cloud one would obtain: 



{n{x))^,^ = {N-i,ip-t\n{x)\N-i,ip-t) 



= N e|m(x;t)r + (l-0|W^2(a:;t)|2 



+ 2^01 - i)ne{w^{x- t)W2{x- 1) 



) 



(17) 



exhibiting the well known interference pattern. More in general, if at t = 0, the time at 
which the confining potential is switched off, the system is in a generic many body state 
described by the density matrix p, the average density after a free evolution of the system 
up to time t, p ^-^ p(t), will be given by 



As always in quantum mechanics, this mean value refers to statistical averages over many 
experimental runs, where each time the system is prepared at t = in the same state p. 
However, as explained before, this is not what it is done in a typical experiment, since the 
direct measure of the observable n{x) is usually problematic. Instead, one takes a picture of 
the released atom cloud at time t, collecting the resulting absorption image, while loosing 
the sample. 

When the two clouds can be considered independent or described by a state like |iV/2; iV/2), 
this procedure is commonly assumed to correspond to a projection onto a fixed-phase state, 
with definite ^ = 1/2 [121 E]- the following, we extend this treatment by considering 
projections of the form |A^; ^, 99; t)(A^; ^, yj; t| onto a state with a definite relative phase 
and amplitude ^, although randomly given, producing the average density (Ti(a;))g^Jl| : our 
purpose is to compute the average over all the absorption images within this approach and 
to compare it with the expression f|T8l) . As it will become clear below, the inclusion also of 
the variable ^ is necessary to preserve the overall probability. Further, if the whole exper- 
iment is repeated, preparing the system in the same state p and taking its picture after a 
free expansion up to the same time t, a different phase and amplitude will be selected. The 
distribution of the obtained values over many repetitions is determined by the initial state p 
through the probability (A^; ^, v?; t|p(t)|-/V; t), which indeed gives the weight in p of the 
state configuration with fixed phase ^ and amplitude ^. 

For large number of atoms A^, this description is perfectly in agreement with what it is 
experimentally measured through the absorption images, that exhibit interference patterns 
irrespectively from the initial state, being either superfluid or Mott. Indeed, in many-body 
physics, one can assimilate ensemble averages with mean values with respect to macroscopi- 
cally occupied many-body states, provided the number of particles involved is large enough 

^This result can be naively understood by interpreting the formation of the absorption image as the result 
of the interaction of the system with a classical, macroscopic measuring apparatus: many atoms concur to 
the formation of a single pixel in the image and this is possible only if all atoms are in a same coherent 
superposition; as explained before, this situation is described by the coherent states ([3|), which in turn are 
the more classical among all quantum states. 
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[21]. Therefore, the larger the number of atoms the system contains, the better a single 
absorption image will model the average of the density operator n{x) in the state |A^; ^, ip; t), 
thus reproducing the interference pattern {n{x))^^^ given in (fT7|) [14]. 

In a typical experiment, all obtained absorption images are subsequently superimposed, 
giving rise to the average 

1 

<^(^)) = Ar5Z , (19) 

i=l 

with A/" the total number of taken pictures. In this sum, any specific density pattern 
(n(x))^.^^., with values C,i and (pi of amplitude and phase, will occur z/j times, i.e. with 
a frequency Pi = Vi/M:, the average f[T^ may then be more conveniently rewritten as 

(n(x)) = 5^ p„ (n(a;))5„,^„ , (20) 

where the sum is now over the set of distinct absorption images. In the limit of large A/", the 
above sum becomes an integral over all possible values of amplitude and phase. Therefore, 
the procedure of averaging over all obtained images really corresponds to summing the mean 
values of n{x) over all states |A^;^,9?;t) with definite ^ and (p, each contribution weighted 
with the associated occurrence probability pa ~ {N] ^, if] t\p{t)\N] ^, ip; t), determined by the 
initial state p. Mathematically, this procedure is then described by the operation of trace of 
n{x) over the density matrix 

m = dip {N;^,ip;t\pit)\N;^,ip;t) \N; ip;t) {N; ip;t\ . (21) 

More precisely, the transformation p \^ p defines a linear map Vat acting on the space of 
all density matrices: it is a realization of a Positive Operator Valued Measure (POVM) 
[TB| [13, dB], a generalization of the more standard von Neumann measurel^ 

p{t) = V^(p(t)) = [ d^ [ d^ V{N; e, t) pit) V{N; e, v?; t) , (22) 







with 



V{N;^,ip;t) = \N;^,^;t){N;^,ip;t\ . (23) 

V Zn 

The proper quantum mechanical expression for the matter density profile resulting from the 
average over all absorption images is then given by: 

{n{x)) = ^— dU dip{N;^,ip;t\pit)\N;^,ip;t) {N;^,ip;t\nix)\N;^,ip;t) , (24) 

^TT Jo Jo 



Strictly speaking, after the release of the trapping potential, the integral Vn+i 

ace 



yj'^^ J^^ d(fV{N;^,(p;t) gives the projector operator on the {N + l)-dimensional subsp. 

spanned by the vectors ^ and not the identity 1 over the whole Hilbert space; therefore, the POVM should 
be more correctly defined by the set {V{N;£^,(p;t), 1 — Vn+i}- However, since all system states at time t 
come from the free evolution of states belonging to this subspace, the action of 1 — Vn+i has no effect, and 
the definition (|22|) given below follows. 
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which is clearly different from the average {n{x))p(t) in (fT8|) . that is usually adopted in 
interpreting the experimental data [31- [S], [12] • 

The difference between these two expressions of the mean matter density may have exper- 
imental relevance. In this respect, the most interesting case is represented by the Fock states 
([2]), describing situations for which the system has a definite number of atoms in the two 
wells: they can be easily realized in practice by raising the barrier between the two traps. 
Let us assume that, just before the release of the trapping potential, the system be prepared 
in a state containing k atoms in the first well and N — k m. the second, so that 

p= |A^;A;)(A^;A;| . (25) 

In the limit of large A^, the distribution of matter in a single absorption image taken after 
a free evolution up to time t is still given by ( fTTI) . with a randomly picked relative phase 
and amplitude ^. Nevertheless, while the distribution of the possible values of (/? over many 
images remains flat, that of ^ follows a time- independent binomial law: 

{N-^,^-t\p{t)\N-^,V-t)= {N-C,r,t\N-k;t) = f Je'(l-0^"'- (26) 

Therefore, if the experiment is repeated many times, each single absorption image will show 
interference fringes with the same spatial period but with randomly distributed offset posi- 
tion, while the amplitude ^ remains essentially constant through all images, since for large 
the distribution fl2Bl) is peaked around its average value. As a consequence, by superim- 
posing all absorption pictures, the interference pattern is averaged away, as it is easily found 
by inserting (l25ll into the expression fl2^ : indeed, by performing the integrals, one explicitly 
finds: 

{n{x)) = j;^[{k + l)\W,{x,t)\' + {N-k + l)\W2{x,t)\'] . (27) 

On the other hand, by adopting the definition (fTSjl . using (l25!l . one would instead obtain 

{n{x))k = k\W^ix,t)\^ + (N -k)\W2ix,t)\^ . (28) 

This two expressions coincide only for equally filled wells, k = N/2; since this is the situation 
encountered in most experimental setups, this explains why in all discussions the definition 
( IT8|l for the average density has been adopted instead of the quantum mechanically correct 
one given in fl24|) . 

Further, notice that the difference between the above two expressions is not vanishing as 
A^ becomes large, 

{n{x)) - {n{x))k ^ {\Wi{x,t)\^ - \W2{x,t)\^) + 0{1/N) , (29) 

and therefore it is of relevance in certain experimental situations. This is the case for states 
with very unbalanced fillings of the two wells, for which k is very small. By selecting a time 
of flight t for which the atom cloud is sufficiently spread out to be visible in the absorption 
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images, but such that the two wave functions Wi{x,t) and W2{x,t) are still sufficiently well 
separated in space, one should be able to experimentally measure the weight in front of the 
two contributions |Vri(a;, t)p and |W2(a;, t)p and thus quantitatively test the validity of the 
assumption of the projection on coherent states during the absorption imaging process. 

As a further remark, notice that if instead the system is prepared at t = in a superfluid 
state, described by the density matrix 



all absorption images taken after a time t will show exactly the same interference pattern, 
described by the given amplitude ^' and phase ip'. This is a consequence of the (large A^) 
orthogonality of the coherent states |A^; ^, ip, ) as given by ([5]), which remains true for all t, 
since the time evolution is unitary. By superimposing all taken images, using and ([5]) 
one then finds 



so that the obtained average density coincides with the mean operator density, thus repro- 
ducing the interference pattern given in (fT7|) . as observed in actual experiments. 

Finally, let us notice that if the double-well trap is ffiled with non-interacting fermions 
instead of bosons, no inteference pattern can emerge in a single absorption image, since 
the fermions can not macroscopically occupy the same state. However, if a small attractive 
interaction is switched on, such that the new ground state is BCS-like with a sufficient 
number of Cooper-pairs, then the description of the system in terms of the Bose-Hubbard 
Hamiltonian ([1]) holds also in this case [251 [261 [211 , ^^^^ ^ consequence, it is expected that 
single absorption images will show inteference fringes. 
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